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Born 2 decades ago in the study of the SUSY breaking mechanism of higher dimensional 
quantum field theories, ^ Supersymmetric Quantum Mechanics (SQM) has so far been 
considered as a new field of research, providing not only a supersymmetric interpretation 
O I of the Schrodinger equation, but interesting answers in all sorts of non-relativistic quantum 
^ I mechanical systems. Particular points to be mentioned include the better understanding 
(-^ • it brought of the exactly solvable, 0-0, the partially solvable, 0], 0, the isospectral, 
D . and the periodic potentials, 0. Recently the association of the variational method 

^ ! with SQM formalism has been introduced to obtain the approximate energy spectra of 
^ I non-exactly solvable potentials, P-[PD| 



> 



Works of reference p| introduce a scheme based in the hierarchy of Hamiltonians; it 



permits the evaluation of excited states for one-dimensional systems. In reference [10] a 
new methodology based in an ansatz for the superpotential which is related to the trial 
CN ! wave function is proposed. This new methodology has been successfully applied to get 
the spectra of 3-dimensional atomic systems and it is illustrated here through systems 



■ well fit by the Hulthen, the Morse and the screened Coulomb potentials, [0, |Tl|], |T2 



As a byproduct of this investigation a new exactly solvable potential has been found, a 
generalization of the Hulthen potential, which presents, in 1-dimension, the property of 
shape invariance. 

Q^' Here, these results in a concise form and a sketch of the new potential are presented. 

(D ■ 

As we learn from the basis of quantum mechanics, (see for instance [|13|)5 the variational 
method was conceived to be an approximative method to evaluate the energy spectra of a 
^ ■ Hamiltonian H and, in particular, its ground state. Its central point is the search for an 
' optimum wave-function \I'(r) that depends of a number of parameters. This is called the 
trial wave-function. The approach consists in varying these parameters in the expression 
for the expectation value of the energy 

until this expectation value reaches its minimum value. This value is an upper limit of the 
energy level. Even though this method is usually applied to get the ground state energy 
only, it can also be applied to get the energy of the excited states. 

Thus most important for the variational method to work is the aquisition of this opti- 
mum wave function. At this crucial point SQM is used to obtain this function. Based in 
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physical arguments, an ansatz for the superpotential is proposed and, through the super- 
algebra, the trial wave function is evaluated. By minimizing the energy expectation value 
with respect to a free parameter introduced by the ansatz the minimum energy is found. 
Consider a system described by a given potential Vi. The associated Hamiltonian Hi 



can be factorized in terms of bosonic operators, in ^ = c = 1 units, fl^- [jT7|. 



Hi = -\^,+Vi{r) = AiAr+E^o^ (2) 

where e'^'^ is the lowest eigenvalue. Notice that the function Vi{r) includes the barrier 
potential term. The bosonic operators are defined in terms of the so called superpotential 
Wi(r), 

As a consequence of the factorization of the Hamiltonian ifi, the Riccati equation must 
be satisfied, 

- W[ = 2Vi (r ) - 2E^^^ . (4) 

Through the superalgebra, the eigenfunction for the lowest state is related to the super- 
potential Wi by 

^W(r) = Nexp{- rWi{f)df). (5) 





What is clear is that if the potential is non-exactly solvable, the Hamiltonian is not 
exactly factorizable, in other words, there is no superpotential that satisfies the Riccati 
equation. The Hamiltonian can however be factorized by an effective potential. In this 
case, the Riccati equation is exact. What we do is to make an ansatz for the super- 
potential using physical arguments that approximate the effective potential to the true 
potential. Having a superpotential we use the superalgebra to evaluate the wave function 
which will definitely depend on a free parameter, the variational parameter. 

We stress that in fact we are dealing with an effective potential Veff that satisfies Riccati 
equation, i.e., 

Y 

where Wi = Wi{fi,) is the superpotential that satisfies (§) for /i = /i, the parameter that 
minimises the energy of eq.(|I|). 

The Hulthen Potential 

The Hulthen Potential, in atomic units, is given by 



Veffiy) = ' ^ ' + Eifi) (6) 



VHir) = -^-z^, (7) 



where 5 is the screening parameter. This potential has been used in several branches of 



Physics, (see [|T8| and references therein). Its Hamiltonian is written as 



^ _ I (P 6e-^'- l{l + 1) 



The associated Schrodinger equation is only solved in closed form for the case / 
is waves) 



0, 



I9| . Other than this case, the potential barrier term prevent us to solve the 
Schrodinger equation and, in SQM superalgebra language, to build the superfamily as can 
be done when / = 0. This case however serves as a basis to construct a superpotential for 
the / 7^ case, [|T0|] , 



Wi{r)=B, 



-Sr 



1-e 



-Sr 



where 



--(l + v/l + 4/(/ + l)). 



2 Bi 



(9) 
(10) 



This in fact defines, through the Riccati equation, an effective potential whose functional 
form is 

, , Se-^^ 1(1 + 1) 5^e-^^' 



1 - e 



+ 



CI 



'ID 



We note that, for small values of 5, the second term of ([TT]) gives us a potential barrier 
term of in first approximation. 

For the state 2p we use equation with / = 1 and evaluate the wave function, changing 
5 by the variational parameter /i, i.e.. 



vl/«(r,/.) 



:i-e 



(12) 



The energy is obtained by minimisation with respect to fi. Thus, the equation to be 
minimised is 

^ _ r^.(r)Hi^-lS^ + ^]^.(0^^ 

and the integration is carried out numerically. Our explicit values for the 2p, {1 = 1) 
energy states for some values of the parameter 6 are listed bellow in Table 1. They are 
shown together with direct numerical integration data. 



State 


Delta 


Variational result 


Numerical Integration 


2p 


0.025 


-0.112760 


-0.1127605 




0.050 


-0.101042 


-0.1010425 




0.075 


-0.089845 


-0.0898478 




0.100 


-0.079170 


-0.0791794 




0.150 


-0.059495 


-0.0594415 




0.200 


-0.041792 


-0.0418860 



Table 1. Energy eigenvalues as a function of the screening parameter for the states 2p, 
[eq.([12|)]. Comparison is made with numerical data of Ref.[|l^]. 



The Morse Potential 



The three dimensional Morse oscillator, suitable to describe a diatomic system, can be 
written as, 

Vm = D{e-^^'^'-''^^ - 2e'<'-^^^) (14) 



where D is the dissociation energy, r^. is the equihbrium internuclear distance and a is the 
range parameter. We rewrite the original Schrodinger equation = E^^ in terms of a 
new variable y, 



J 



(15) 



where y = ar and the constants are set like 



ye 



2mD 



2^2 



E 



h a 



and the parameter m is the reduced mass of the molecule. 



(16) 



For the case / = the system is again exactly solvable and it is used to provide 
information to the ansatz to be made to the superpotential for the I ^ case, for which 
an analytical exact solution cannot be determined. The superpotential takes the form 



Wi(y) = -\e-^y~y^^ - ^^-^ + c. 

y 



(17) 



The argument is that the first term is taken from the one- dimensional resutls, ||21I], the case 
of / = 0. The knowledge of the second term comes from the study of three-dimensional 
potentials, |2^. The related effective potential is given by 



eff 



/+1 



/+1 



Em) 



The eigenfunction obtained from eq. ([TTp is then 



\E'(?/) oc e 



.Xe-(y-ye) yl+l ^-Cy 



(19) 



Using this expression as a trial wavefunction in the variational method we change the 
parameter C by the variational parameter n, i.e.. 



*(?/,C = /i) oc e 



(20) 



The energy is then obtained by minimisation of the energy expectation value with respect 
to 11. Thus, the equation to be minimised is 



E„ 



dy 



^f.{y)dy 



(21) 



!^^,{yYdy 

We have used this expression to minimize the energy expectation value of various 
molecules: H2, HCl, CO and LiH, [|lT|. The explicit values of the energy for n = 
and different values of I are shown bellow for the H2 molecule, for known values of their 
respective potential parameters, p3l: D, a, and m. 



1 




Shifted 1/A^ 
pxDririsiori rpsiilts 


Modified shifted 1/A^ 
pxnansion rpRiilts 

V_j l-J CAj X X kj7 X XX X V_j kjj I. -I. X U 


T^;ic/^pt INFnrnprip^il 

J /vVOj U J. 1 Li-XXXVjX X^.^CtjX 





-4.4758 


-4.4749 


-4.4760 


-4.4759 


5 


-4.2563 


-4.2589 


-4.2592 


-4.2589 


10 


-3.7187 


-3.7247 


-3.7252 


-3.7242 


15 


-2.9578 


-2.9663 


-2.9670 


-2.9654 


20 


-2.0735 


-2.0839 


-2.0846 


-2.0826 



Table 2. Energy eigenvalues (in eV) for different values of / for H2 molecule, with 
D = 4.7446eV^, a = 1.9426A~^, = 0.7416A and m = 0.50391amM. Comparison is made 
with results from ref. . 

The Screened Coulomb Potential 

The screened Coulomb potential is given, in atomic units, by 

„—Sr 



sc 



(22) 



where 6 is the screened length. The associated radial Schrodinger equation includes the 
potential barrier term and it is given by 



1 rf2 e-^' i{i + iy 



+ 



2r2 



(23) 



where the unit length is h'^ /me^ and the energy unit is eg = —me^/Ti^. 

In order to determine an effective potential similar to the potential in the Hamiltonian 
(P3|), that is the screened Coulomb potential plus the potential barrier term, the following 
ansatz to the superpotential is suggested 



5e 



-5r 



1 



(24) 



(/ + 1) 2 
Substituting it into (H), one gets 

^o(r) = (1 - e-^'-y+^e-^^-^^'. (25) 

Assuming that the radial trial wave function is given by (P^D, replacing 6 by the variational 
parameter /x, i.e., 

^^(r) = (1 - e-^^)'+ie"((^"^)'-^ (26) 
the variational energy is given by 



(27) 



Thus, minimizing this energy with respect to the variational parameter /i one obtains the 
best estimate for the energy of the screened Coulomb potential, |1^ . 



As our potential is not exactly solvable, the superpotential given by eg. does not 
satisfy the Riccati equation but it does satisfy it for an effective potential instead, 

Veff 

VefAr) = ' ^ ' + Em) (28) 
where Wi = Wi{5 = Ji) is given by eq.(^) and p, is the parameter that minimises the 



energy expectation value, p7|). It is given by 



be- 



ef f\ 



+ 



/(/ + 1) 



-2&T 



1 



1 



2' 



(29) 



l-e-^"- ' 2 (l-e-^^)2 2'/ + l 
where S = ft that minimises energy expectation value. One observes that for small values 



of S the first term is similar to the potential (^2]) and the last is approximately the po- 
tential barrier term. This observation allows us to conclude that the superpotential (p^ ) 
can be used to analyse the three dimensional screened Coulomb potential variationally 
through the trial wavefunction (pSj). 





2] 


□ 


3d 


4f 


6 


Variational 


Numerical 


Variational 


Numerical 


Variational 


Numerical 


0.001 


-0.2480 


-0.2480 


-0.10910 


-0.10910 


-0.06051 


-0.06052 


0.005 










-0.52930 


-0.05294 


0.010 


-0.2305 


-0.2305 


-0.09212 


-0.09212 


-0.04419 


-0.04420 


0.020 


-0.2119 


-0.2119 


-0.07503 


-0.07503 


-0.02897 


-0.02898 


0.025 


-0.2030 


-0.2030 


-0.06714 


-0.06715 






0.050 


-0.1615 


-0.1615 


-0.03374 


-0.03383 






0.100 


-0.09289 


-0.09307 











Table 3. Energy eigenvalues as function of the screening parameters S for 2p (1 = 1), 
3d {I = 2) and 4/ (/ = 3) states, in rydberg units of energy. Comparison is made with 
results from references |P4|| -||29||. 

Comments 

We have proposed trial wavefunctions to be used in the variational calculation in order 
to determine the energy eigenvalues of some atomic systems described by the Hulthen, 
the Morse and the screened Coulomb potentials. These functions were induced from the 
formalism of SQM. Using physical arguments it is possible to make an ansatz in the 
superpotential which satisfies the Riccati equation by an effective potential. The superal- 
gebra enables us to take this superpotential and to evaluate the trial wavefunctions that 
contains the variational parameter, the parameter that minimises the energy expectation 
value of the energy. The results are in very good agreement with the ones available in the 
literature. 

Thus this new methodology that associates the variational method with SQM is a sim- 
ple and good alternative procedure that enables the evaluation of the energy spectra with 
reasonable accuracy. In particular, when applying the approach to the Hulthen potential 
it was found that its effective potential was linked to a new exactly solvable potential, 
that presents, in one-dimension the property of shape invariance. At that time the focus 



of our attention was not on the properties of this new potential, whose main results are 
presented bellow. 

A new shape-invariant potential 

As mentioned above, in [|l^] when dealing with the evaluation of the trial wave function, 
the following one-dimensional potential, a generalization of the Hulthen potential, was 
found. 



n^) = ^ .. .-s.,2 -Br—^- (30) 



Notice that when setting A = and B = 1 the original Hulthen potential is recovered, 
eq. d^). For convenience we redefine the constants A and B such that V{x) becomes 

Vi{x) = a,{a, - ^) 2(i_g-fa)2 ' «i(2^i + 2(1 - e"'^^') " ^^^^ 
From the basis of SQM the following superpotential 

W, = -a,^—-^ + h (32) 
factorizes the related Hamiltonian 

Hi = -\^, + ^i(^) = + E'i^ (33) 

and satifies the Riccati equation associated 



Wl{r)-W[{x) = 2Vi-2E^^'^ 



where the lowest energy-eigenvalue is given by 

4" = -f (34) 
From the superalgebra, the associated eigenfunction, is given by 

^« = (l-e-^")^e-''i". (35) 

The condition we must impose on the above wave-functions is to vanish at infinity and at 
the origin, i.e., 

^>0, 6i>0 ^ ai,bi,6>0 (36) 



This process of factorization can be repeated and the whole hierarchy of this potential 
can be evaluated. The general form for the superpotential is written as 

g-5x 

Wn+l = -Qn+l ^ _ + K.+l (37) 



which are related to the potentials 



with lowest levels given by 

and with the constants being given by 

ttn+i = ai + nS (40) 

and 

= 7r^(ai(26i + 5) - 25(nai + _ ^. (41) 

For the particular case where n = and fixing the constants to 

ai = 5 , 26i + 5 = 2 (42) 
the original Hulthen potential is recovered, 



which can also be factorized, as already known, ||T0|. Its whole hierarchy of potentials is 
given by 

_n{n- l)52e-25^ 5(2 + n{\ - n)5)e-^^ 
"~ 2(l-e-^-)2 2(1 - e-^-) ^ ^ 



with lowest energy levels 
and related superpotentials 



4"'^-i(i-^y^ (45) 



Wn{x) = -flnj^^ + bn (46) 

where the lowest states are 

^J,") = (1 - e-5^)an/'5e-fena; (47) 

with constants given by 

ttn/S = n , bn = — . (48) 

n 2 

At this point we remark that the Hulthen potential is non-shape- invariant, [0]. This is 
clear by observing the n = 1 and n = 2 cases of equation @), whereas the new potential 
clearly is shape invariant, since all the potentials preserve the shape in the hierarchy eq. 
(^). More details shall soon be found in . 
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